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ABSTRACT 

The evolution of the abundance of galaxy clusters is not a reliable measure of f2 if there 
are features on scales of a few Mpc in the primordial power spectrum. Conversely, if 
we know the cosmological model parameters from other measurements, the cluster 
abundance evolution permits us to probe features in the power spectrum that are in 
the nonlinear regime at the present epoch, and hence difficult to discern directly from 
current epoch measurements. 

We have investigated the influence of an artificially introduced Gaussian fea- 
ture on an otherwise unperturbed SCDM power spectrum on scales corresponding 
to k ~ 0.4 — 0.8 /jMpc . Using these modified spectra as an input to cosmological 
./V-body simulations, we are able to show that in terms of the cluster abundance evo- 
lution, a SCDM model displays characteristics similar to an OCDM model. However, 
strong modifications would also be visible at a redshift z = in the dark matter power 
spectrum whereas minor alterations to the usual SCDM spectrum are washed away 
by non-linear evolution effects. We show that alterations to the dark matter power 
spectrum like those presented in this paper do not leave any imprint in the present 
density fluctuation spectrum and the velocity distribution of galaxy clusters; nearly 
all models agree with each other and do not coincide with our fiducial OCDM model, 
respectively. We therefore conclude that features with characteristics such as discussed 
here might not be detectable using observations of the galaxy power spectrum, the 
local cluster abundance or the large-scale velocity field as measured by the velocity 
distribution of galaxy clusters. 

The only quantity that shows a pronounced difference at the present epoch be- 
tween our models under investigation is the halo-halo correlation function which ap- 
pears to be strongly biased with respect to an unmodified SCDM model. This is due to 
a lack of power on certain scales which subsequently modifies the relative amplitude of 
high- and low- A; waves. Apart from observations of the evolution of cluster abundance, 
measurements of the Lyman a forest at high redshift could put constraints on possible 
features in the power spectrum, too. 

Key words: large scale structure - cosmology: theory - cosmology: large scale struc- 
ture of Universe 



1 INTRODUCTION 

The distribution of matter on large scales in the Universe 
is supposed to have evolved by gravitational interactions 
from seeds originating in quantum fluctuations that were 
stretched to cosmological dimensions by inflation. Ordi- 
narily, this initial spectrum of fluctuations is assumed to 
be a featureless power law. This requires the least num- 
ber of parameters and is produced in the simplest models 
of inflation. However, different observations, such as clus- 
ter redshift surveys (Einasto et al. 1997) and galaxy sur- 
veys (Hamilton & Tegmark 2000, Gaztanaga & Baugh 1998, 
Broadhurst et al. 1990) reveal possible traces of features in 



the matter power spectrum. In addition, a lower than ex- 
pected second peak in the CMB spectrum as measured by 
BOOMERANG (Lange et al. 2001) appears to be incom- 
patible with the simplest models of cold dark matter and a 
scale invariant primordial spectrum. 

While the results, particularly from the surveys are still 
not beyond statistical doubt, a genuine feature(s) in the pri- 
mordial spectrum cannot be ruled out. Several mechanisms 
have been proposed that could generate features in the pri- 
mordial spectrum during the epoch of inflation. They com- 
monly involve an extension of the simplest one field infla- 
tion model, e.g. by coupling the inflaton to a massive par- 
ticle (Chung et al. 1999) or considering two-field inflation 
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Figure 1. CDM input spectra at redshift 2 = 0. 

(Lesgourgues, Polarski & Starobinsky 1998). There are also 
more exotic ideas (see e.g. Martin, Riazuelo & Sakellariadou 
1999). Another idea proposed is to Taylor expand the pri- 
mordial power spectrum to include higher order terms that 
account for a running spectral index (Lidsey et al. 1997, 
Hannestad, Hansen & Villante 2000). The latter is able to 
introduce a very broad negative or positive bend into the 
power spectrum. 

While possibly providing a motivation for observed fea- 
tures on a specific scale (i.e. scales of about 100 Mpc), the 
inflationary mechanisms proposed appear in principle capa- 
ble of producing features on other scales as well. Here we 
examine the effect of primordial features on smaller scales 
k ~ 0.4 — 0.8 /tMpc -1 , where due to (the onset of) non- 
linear evolution and the problem of biasing, a connection 
with the primordial spectrum is much harder to establish. 
We have previously examined the effects of such bumpy 
power spectra on the cosmic microwave background (Grif- 
fiths, Silk and Zaroubi 2001). 

In Section ^| we will present the modifications applied 
to an otherwise unperturbed SCDM power spectrum. These 
spectra were then used as input to cosmological iV-body 
simulations as described in Section ^[ Section ^ deals with 
the complete analysis of these simulations with respect to 
the evolution of the power spectrum, velocity statistics, the 
masses of galaxy clusters, the halo-halo correlation function, 
and density profiles for a selection of halos. In Section H 
we try to link our numerical results to analytical predic- 
tion mainly based on the Press-Schechter theory (Press & 
Schechter 1974). We close with a discussion of our main re- 
sults in Section H. 



2 THE POWER SPECTRA 

We focus our attention on an SCDM model with parameters 
as given m Table |l} Although in this paper we are looking 
at properties of clusters and particularly their abundance 
evolution, we decided to use only COBE-normalised spectra 
with spectral index n = 1 for modifications. This is because 
the features we add to the spectra are on scales that at the 
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Figure 2. Evolution of CDM spectra. 



present time are subject to non-linear evolution and thus 
are expected to lead to a deviation from the erg — Q, relation 
(Eke et al. 1996) that would otherwise be used to cluster- 
normalise the spectra. We also focused our attention mostly 
on negative amplitude features - or dips - in the spectrum 
of SCDM models. In addition to the significance of features 
for structure evolution on relatively small scales, dips will 
reduce the overall power in the spectrum, pushing the nor- 
malisation and possibly other properties towards what is 
expected in OCDM models. 

The corresponding SCDM power spectrum was calcu- 
lated using the publicly available CMBFAST code (Seljak 
& Zaldarriaga 1996) and hence is COBE-normalised. For 
Model lb we lowered the amplitude to reach a normalisa- 
tion that agrees with the cluster abundance as described in 
Eke et al. (1996). All other models are based on Model la 
with no further modification to the amplitude than intro- 
duced by the artificially added feature. The modified power 
spectra follow the equation: 

P mod (fc) = P(k) ■ (1 ± Aexp[-0.4( logfc ~ logfco ) 2 ]) (1) 

<^mod 

where P(k) is the unmodified spectrum. 

These features are completely Gaussian in log-space and 
can be described by their width <r mo d, height A, and the lo- 
cation fco of the bump (+) and dip (-), respectively. Their 
parameters are also summarised in Table |l|. Each modifi- 
cation (characterised by the Model no. 2-6) was applied to 
the SCDM spectrum corresponding to Model la keeping the 
COBE normalisation fixed. A visual impression of the result- 
ing power spectra can be found in Fig. |l] where all spectra 
are plotted (linearly extrapolated to redshift z = 0.0). These 
spectra are now used as an input to our initial conditions 
generator for the cosmological iV-body simulations. 

Again, we decided to use a COBE-normalised SCDM 
power spectrum as reference model because our intention 
was to investigate the influence of artificially introduced 
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Table 1. Specifications of the artificial modifications to the un- 
derlying Gaussian power spectra. The added features are of log- 
normal form, Afe^, a 2 ]. All modifications were superimposed onto 
an otherwise unmodified, COBE normalised SCDM model (Model 



la). The corresponding mass scale M 
given in the last column. 
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features on scales corresponding to galaxy clusters (i.e. 
8/t" 1 Mpc). For this reason we chose the COBE normali- 
sation even though this SCDM model seems to be rather 
unattractive or even ruled out nowadays. However, by tak- 
ing away power on cluster scales (as done for the majority of 
our 'feature' models) the value of as drops as we do not ap- 
ply any re-normalisation of the power spectra. The purpose 
was not to find a new standard model which fits the ob- 
servational data better but rather to analyse the influence 
of such features on an interesting range of cluster quanti- 
ties; and this is more easily followed in a more distinctive 
structure formation scenario. 



3 THE TV-BODY SIMULATIONS 

The simulations were carried out using Couchman's AP 3 M 
code (Couchman 1991). All simulations were performed with 
128 3 particles in a box of side length 256ft. -1 Mpc; the (co- 
moving) force resolution was fixed at 100/i -1 kpc for all 
runs. We evolved the particle distribution from redshift 
z = 30.0 until z — 0.0 in 5000 steps. The box size was 
chosen such that the scales on which the Gaussian features 
had been added lie comfortably within the range covered 
by the simulations, e.g. our simulations cover the k-range 
from fc m i n = 0.0245 (limit set by box size) to fc ma x = 1.571 
(limit set by particle number for representing the waves ini- 
tially present), and the modifications lie clearly within that 
range (cf. Fig. |^). The particle-particle summation part of 
the AP 3 M code guarantees that we properly follow the evo- 
lution of all initially present waves especially for models 2-4 
(cf. Couchman 1991). 

For identifying particle groups within our numerical 
simulations we used the standard friends-of-friends algo- 
rithm (Davis et al. 1985) with the linking-lengths 11 = 0.2 for 
the SCDM, and 11 = 0.17 for the OCDM model, respectively 
(cf. Knebe & Miiller 1999). 



4 ANALYSIS 

4.1 Power Spectrum Evolution 

In Fig. |^ we plot the evolution of the dark matter power 
spectrum for models 1 through 6 (Model lb and the OCDM 
model are left out for clarity as there is nothing unexpected 
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Figure 3. Velocity distribution for all dark matter particles 



to observe). We see that the spectrum appears to be con- 
verging towards the unperturbed spectrum at late times. 
The feature nearly vanishes completely, leaving no further 
imprint in the power spectrum. The thin solid line cross- 
ing the whole plot corresponds to the unmodified, COBE- 
normalised SCDM spectrum linearly extrapolated to z = 
0.0. Comparing it with the spectra derived from the numer- 
ical simulations, we notice a significant boost in power on 
small scales due to non-linear evolution. However, Model 6 
with the very prominent dip at 16 Mpc nearly matches 
the linearly extrapolated SCDM power spectrum: at late 
times the non-linear boost of power almost exactly compen- 
sates the lacking power at the location of the dip. 

We may conclude that even when starting with a promi- 
nent feature on small scales, there might only be little if 
any evidence for it left in the present universe when looking 
at the dark matter power spectrum P(k). In any case, the 
evolution of the models has to be different, and one needs 
to think of other ways to investigate the influence of such 
features in the present day universe or to detect them in 
the evolution of P(k) and related quantities. As far as the 
spectrum itself is concerned current galaxy clustering sur- 
veys offer little prospect of detecting these features in the 
evolution of P(k). Observations of the Lycv forest, however, 
appear to be a promising tool for constraining P(k) at high 
redshift (see e.g. Croft et al. 2000, Weinberg et al. 1998), 
particularly on the scales under investigation here. 
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Figure 4. Velocity distribution for particle groups identified us- 
ing a standard friends-of-friends algorithm with linking length 
0.2. Only particles groups heavier than 2 • lO 13 /!" 1 Mq are taken 
into account. 



4.2 Velocity Statistics 

Next we consider the distribution of velocities for a) the dark 
matter particles and b) friends-of-friends particle groups. 
The peculiar velocity field results from the gravitational ac- 
celeration that develops from initial density fluctuations in 
the early universe. Clusters of galaxies can therefore be used 
as tracers of the large-scale peculiar velocity field (Bahcall, 
Cen & Gramann 1994). 

In Fig. ^| we show the probability distribution P(v) of 
dark matter particles with peculiar velocities in the range 
v ± dv along with the integrated distribution P(> v). The 
same is plotted for friends-of-friends groups identified in all 
runs in Fig. ^| All curves are normalised by the total number 
of particles and total number of particle groups in the respec- 
tive model. Here we also show the data for SCDM Model lb 
as well as the fiducial OCDM model. Even though there 
are differences of a factor of two at the high velocity end of 
the distribution, it might be difficult to observationally dis- 
tinguish these differences between modified and unmodified 
models. 

The deviations between the cluster-normalised SCDM 
model and the COBE-normalised one are bigger than any 
differences between our 'feature' models. This deviation in 
velocities tending to lower values in the cluster-normalised 



model can be explained by the lower value for erg: the am- 
plitudes of the initial density fluctuations are smaller and 
hence it takes longer to accelerate particles (and clusters) 
to high velocities; in the course of a simulation the veloc- 
ity distribution function is always similar to a Maxwellian 
distribution whose peak gradually moves from low to high 
velocities. Apart from this we get increasingly more high ve- 
locity particles leading to a bigger 'tail' in the distribution. 
However, the most prominent difference is in the high veloc- 
ity tail of the distribution for clusters. Only here it is possible 
to discriminate between models as our dip models show a 
significant drop in the integrated probability. We do not ob- 
serve any galaxy clusters with peculiar velocities higher than 
3000 km s _1 whereas there is a distinctive number of these 
objects in the unmodified Model la (and the bump Model 2). 
A detailed check showed that this tail is mainly due to ob- 
jects with masses M <J 10 14 /i _1 M©. When only taking into 
account particle groups with heavier than 10 14 /i _1 Moall 
curves for the 'feature' models fall on top of each other. 
However, the location and the width of the feature does not 
seem to have any significant influence. 

4.3 The Masses of Galaxy Clusters 

The most basic property of a galaxy cluster is its mass M. 
Nevertheless, this quantity can provide a lot of information 
especially when using the (cumulative) distribution of ob- 
jects with a certain mass M; and the evolution of the abun- 
dance of massive clusters within a given mass range is indeed 
one of the corner stones of the currently favoured ACDM 
model (Bahcall et al. 1999). 

4-3.1 Mass Function and Press- Schechter Prediction 

In Fig. [] we show the cumulative mass function of par- 
ticle groups for all our models at a redshift of z = 0.0. 
Even though we could not find well-pronounced imprints 
of the modifications in the power spectra at redshift z = 0.0 
(cf. Fig. [|), we clearly see differences in the amplitude and 
slope of the cumulative mass function n(> M). These de- 
viations are mainly at the low mass end of the resolvable 
mass range, where a positive feature leads to an excessive 
number of objects and a negative feature to a lack of groups. 
This is in general agreement with an excess/lack of objects 
on mass scales corresponding to where the features are lo- 
cated. However, we always end up with the same number 
of galaxy clusters for masses M > lO 15 ^ 1 M (besides for 
Model lb which shows too few massive groups due to the 
low normalisation) . Unfortunately we are not able to resolve 
particle groups lighter than about 2- 10 13 fe _1 M© and more- 
over, the simulation volume is still too small to get a sta- 
tistically significant number of objects heavier than about 
3- 10 15 /i _1 M©. 

To check the validity of our results and get an im- 
pression of how they might generalise for larger/smaller 
mass objects we performed Press-Schechter (PS) calcula- 
tions (Press & Schechter 1974) of the abundance of gravita- 
tionally bound objects. The (differential) number of objects 
for a given mass M can be calculated using their formula: 

dn r2 p 6 C d\no M ( &l \dM 

dM dM= ]JnM^ dhTM 6XP TT (2) 
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Figure 7. Evolution of SCDM cluster abundance. 
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Figure 5. CDM mass function for redshifts z = 0.0 from simu- 
lations (histograms) and Press-Schechter prediction (thin lines). 




with the variance o~m defined as follows: 

om = J P(k)W 2 (kR)k 2 dk , (3) 

where W(kR) is the window function (top-hat in our case) 
for filtering fluctuations in the power spectrum on scales 
characterised by R and hence mass M = 4nR 3 /3. Using 
our model power spectra from Table |l| (plotted in Fig. [jj) 
together with Eq. (g) and Eq. (g) we are able to compare our 
numerically achieved mass functions with the PS prediction. 
The results are shown in Fig. |H| as thin solid lines. 

We observe a similar phenomenon as already seen in 
other comparisons of PS-predicted and iV-body mass func- 
tions (Efstathiou et al. 1988, White, Efstathiou & Frenk 
1993, Gross et al. 1998, Governato et al. 1999, Jenk- 
ins et al. 2001): the PS theory tends to show too many low 
mass objects (about a factor of 1.2) and too few high mass 
objects (again a factor of about 1.2). Apart from that the 
simulations agree fairly well with the Press-Schechter predic- 
tion, even in the cases where we modified the power spectra. 

Anyway, to allow for better comparison between the 
individual models and the effect of the features on n(> M) 
we also plot all mass functions derived from the numerical 
simulations (as already presented in Fig. ^|) in one single 
Figure ^| 

4-3.2 Evolution of the Cluster Abundance 

The evolution of the cluster abundance has proven to be (po- 
tentially) one of the key constraints on the density parameter 
fto (e.g. Eke et al. 1996, Bahcall et al. 1997, Eke et al. 1998, 
Bahcall et al. 1999). This leads immediately to the question 
of how our modifications to the dark matter power spec- 
trum influence this important issue. We have already seen 
that while we might not find hints in the observational power 
spectrum for the features under investigation, there are dif- 
ferences in the mass functions n(> M). In Fig. ^ we show 
the evolution of galaxy clusters with mass M greater than 
lO 14 ^ 1 M . 

The first thing that catches the eye is the difference be- 
tween the COBE and cluster-normalised SCDM model when 
compared to the OCDM model. Only the cluster-normalised 
SCDM model shows a very steep evolution of the cluster 
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redshift z 

Figure 8. Evolution of SCDM cluster abundance normalised to 
unity at redshift 2 = 0. 

abundance whereas the evolution in the COBE-normalised 
one agrees more or less with the OCDM model (cf. Bahcall, 
Fan & Cen 1997, Fig. 2 in their paper). It is now interesting 
to check whether our features affect only the normalisation 
of these evolutionary curves or also the slope. As we can see, 
the variations in the amplitude of the dip (or bump) mainly 
shift the cluster evolution up and down, and the location 
has an influence primarily on the slope of the curve. 

This can be understood in the following way. Features 
on the 8/i -1 Mpc scale are 'dynamically' important, in that 
this scale has evolved from the linear regime (z ~ 1) into 
quasi-linear or non-linear regime by now. Since the rela- 
tive abundance of objects increases dramatically in the non- 
linear regime, we consequently expect still more (less) ob- 
jects at late times if power had been added (subtracted) on 
the corresponding scale, i.e. features on this scale are ex- 
pected to affect the shape of the abundance evolution of 
corresponding mass objects. Conversely, features on scales 
that up to present have remained in the linear regime only 
add (subtract) power on dynamically unimportant scales, 
thus only affecting the normalisation. Q 

Moving the dip from the non-linear scale of 8/i _1 Mpc 
to the semi-linear 16/i _1 Mpc mainly results in an overall 
lack of power and therefore we expect the cluster evolution 
to be closer to the cluster-normalised (erg = 0.52) Model lb. 
However, the evolution for Model 4 agrees fairly well with 
the OCDM Model. This can be seen even better in Fig. ^ 
where we normalised the cluster abundance at redshift z = 
to unity for all models to allow for better comparison of the 
slopes of the curves. 

A similar effect was already observed by Bar- 
riga et al. (2000) where the influence of a step-like feature in 
the primordial power spectrum was studied in the context 
of phase transitions during inflation. However, their features 
are on larger scales and thus would leave traces in current 
epoch measurements of the CMB and large scale structure. 



* Though, if the difference in power is very large (cf. abundance 
evolution for erg = 1.18 and erg = 0.52) there is a significant 
change in the slope of the abundance evolution, too. 
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Figure 9. Evolution of erg in all models. The evolution in all 
models is normalised to be unity at redshift z = 0. 

4.4 erg Evolution 

As the number of objects formed depends on the normal- 
isation of the input power spectrum measured via erg (cf. 
Eq. ^|), it is interesting to check the evolution of this quan- 
tity with redshift, too. Fig. ^ shows this evolution for all 
models. Again, all curves are normalised to the value at red- 
shift a = to allow for better comparison of the slopes. We 
observe no overlap of any feature model with the OCDM 
model: regarding the evolution of the dark matter fluctua- 
tions measured in 8h~ x Mpc spheres (rather than the cluster 
evolution) the OCDM model still differs significantly from 
all other models. Moreover, it is difficult again to distinguish 
between our feature models as well as the differently nor- 
malised SCDM models itself. Only Model 6 shows a steeper 
evolution from redshift z ~ 1 to z ~ 0.5. Hence erg-evolution 
does not provide a discriminant of features in the power 
spectrum. However, Robinson, Gawiser and Silk (2000) have 
shown that erg and cluster abundance combine to probe non- 
gaussianity. 

4.5 Halo-Halo Correlation 

In Fig. [To] we present the halo-halo correlation function 
for our models. We fixed the number density to n = 
5 ■ 10 -4 /i 3 Mpc -3 which actually means to only use the 
N — n ■ V most massive halos with V = 256/i~ 1 Mpc being 
our simulation volume. 

We now observe more obvious differences between these 
models, which can be fully ascribed to our artificial modi- 
fications of the initial power spectrum. The amplitude of 
the cluster correlation function is actually insensitive to the 
amplitude of fluctuations in the density field (Croft & Ef- 
stathiou 1994), which is clearly reflected when comparing 
Model la and lb. Moreover, there are also only moderate 
changes in the cluster correlation function when varying the 
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Figure 10. Halo-Halo correlation function for particle groups. 
The number density of objects was fixed in all models to n = 
5 ■ lCT 4 h 3 Mpc- 3 . 

cosmological parameters from our SCDM Model la/b to the 
OCDM model (cf. Martel & Matzner 2000). But as we move 
on to our feature models not only the amplitude of the cor- 
relation function increases (for the dip models), but also the 
slope changes (i.e. Model 6). This suggests a strongly (and 
non-linearly) biased formation of galaxy clusters in those 
models compared to a 'normal' SCDM model. 

As there was no obvious biasing observed in the dark 
matter power spectrum (which is nothing more than the 
Fourier transform of the dark matter correlation function), 
the effect seen in Fig. [u] is entirely due to a different cluster 
formation scenario in our dip models. This is caused by the 
fact that by adding (subtracting) power only on a certain 
scale the relative strength of different waves becomes more 
and more important. If we put less power into high-fc waves 
(i.e. Model 3 and 4) but leave the largest waves unmodified 
we eventually bias the formation process of galaxy clusters 
with respect to an unmodified SCDM model as can be ob- 
served in Fig. Massive clusters indeed are highly biased, 
and the extent to which one may have difficulty in account- 
ing for this in standard low density models may constitute 
the strongest signature of a possible feature at redshift z — 
in a SCDM model. 

4.6 Density Profiles 

We have already seen that we introduced a strong bias in 
the cluster formation process by altering the ratio of high- to 
low-fc waves amplitudes. This immediately raises the ques- 
tion of whether this will subsequently lead to deviations 
in the shapes of the clusters themselves. We have there- 
fore calculated the density profiles for massive FOF groups 
(M > 10 15 /i _1 Mq) and show the results for a representa- 
tive high-mass halo M ~ 3.5 • 10 15 /i _1 Mq (~ 1700 parti- 
cles) in Fig. [ll]. The curves for models 2 through 6 have 
been shifted upwards by successive factors of 2, whereas the 
OCDM profile was lowered by a factor of 2. 

We can see that modifying the power spectrum does 
not affect the slope and amplitude of the density profiles. 




100 1000 
r [h~~ kpc] 

Figure 11. Density profile for one halo in all models. 

Only the SCDM model lb deviates from the other curves as 
the corresponding halo consists of only about 500 particles. 
Due to the resolution limits of our simulations, we are not 
able to quantify the substructure content of these galaxy 
clusters in more detail as the most massive objects contain 
'only' 2000-2500 particles. But a visual comparison of the 
particle distributions in those clusters provided no obvious 
differences. 



5 ANALYTICAL ESTIMATES 

We have seen above that SCDM models with features at 
small scales may display an evolution of their cluster abun- 
dance that is similar to that in the OCDM model to within 
some constant offset in amplitude. The question remains of 
how far we can push this agreement and to what extent this 
also holds for cluster masses larger than the ones we were 
able to investigate numerically. To this end we performed 
Press-Schechter (Press & Schechter 1974) calculations to 
investigate a larger set of models. However, since we are 
considering the evolution of cluster abundances, this also 
requires a relation for the redshift dependence of the criti- 
cal overdensity S c , which is left unspecified in the original 
PS approach. We determined approximate relations for the 
unmodified SCDM and OCDM models by matching their 
cluster abundance evolution to that of the corresponding 
N-body simulations: 

S^ CDM (a) = 1.2 + 0.6a, 5° CDM {a) = 1.45 + 0.3a (4) 

where a — 1/(1 + z) is the cosmic scale factor at redshift z. 

For the OCDM model the S c relation almost exactly 
reproduces the abundance evolution. The corresponding re- 
lation for the SCDM model was rather chosen to fit both 
unmodified SCDM and featured SCDM models reasonably 
well, matching the N-body evolution to within 10 percent for 
the unmodified SCDM model, typically to within 20 percent 
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Figure 12. Comparison of evolution of clusters of mass M > 
lO 14 M0fe — 1 in simulations and PS using the relations for <5 C from 
the simulations. 



Figure 14. PS calculations of abundance evolution for OCDM 
and SCDM with a broad dip at 10 Mpc h~ 1 . For the mass cuts 
the two agree to within less than a factor of 2 and 3 respectively 
across the whole redshift range. 
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Figure 13. Comparison of abundance evolution for clusters with 
mass M > 10 14 MqIi~ 1 in simulations and PS. All PS calculations 
use S c (a) = 1.2 + 0.6a (eq (§)). 



for models 2 to 5 and 30 percent for model 6. This is shown 
in figures [l^ and [l3| An accurate fit to the unmodified 
SCDM model only would have required 5 c (a) = 1.25 + 0.5a. 

However, the agreement between N-body and PS data 
should be viewed with care, particularly when large sharp 
features are included: At the location of the features, we es- 
sentially introduce a more rapidly changing spectral index. 
This might contribute to a larger discrepancy between N- 



body and PS data. To our knowledge agreement between 
PS and numerical simulations has so far only been estab- 
lished for power spectra with constant (power law) or slowly 
varying (e.g. SCDM) spectral indices - it is not obvious why 
this should also hold when a sharply varying index is being 
introduced, particularly on scales that only now turn non- 
linear. In the following, however, we will assume that this 
is not a problem and use the <5f CDM relation from equation 
(^) as a best estimate for PS predictions of the abundance 
evolution for other small scale features in the SCDM model. 
We also assume that the <5 C relation for OCDM describes the 
abundance evolution in higher mass cuts as well as it does 
for M > 10 14 h' 1 M in our simulations. 

As seen above, one might expect that adding (sub- 
tracting) power only at specific scales leads to a larger 
(smaller) abundance of objects with mass corresponding to 
these scales. This is not obvious for scales that are already 
in the non-linear regime, as we have a coupling across scales, 
and also objects not only get newly created but now have 
also been partly incorporated into larger objects, i.e. the ef- 
fect of a feature on a specific scale spreads out, not only in 
the evolution of the power spectrum as we have seen above, 
but also in the abundance of objects of corresponding mass. 
However, for broad features we may assume that the latter 
is less important, since a whole range of scales and corre- 
sponding masses are affected. 

To check this, we placed a broad Gaussian dip in the 
SCDM model (A = -0.55, cr mod = 2.5) at 10/i _1 Mpc, 
which is expected to affect a similarly broad range of masses 
centred on M ~ 10 14 - 1O 15 M /i _1 . A broad feature of this 
kind is naturally generated by e.g. including a running spec- 
tral index term when expanding the primordial perturbation 
spectrum (Hannestad, Hansen & Villante 2000). 



© 0000 RAS, MNRAS 000, 000-000 



Bumpy Power Spectra and Galaxy Clusters 9 



Comparing Fig. [12| with Fig. [14j confirms that the 
abundance of masses in these ranges (M > 10 14 , M > 
8-10 14 /i _1 M ,the mass range covered by e.g. Bahcall & Fan 
1998) is significantly suppressed. What is more, the abun- 
dance evolution for the mass cuts agrees in logarithmic slope 
and magnitude with that in the OCDM model to within a 
factor of less than two and three across all redshifts up to z 
= 1 for the lower and upper mass cuts respectively. 

Given that observations of the most distant massive 
clusters currently only constrain the abundance to within 
one or two orders of magnitude (see e.g. Bahcall & Fan 
1998), SCDM with broad features towards small scales, 
mimicking (in shape and absolute magnitude) the abun- 
dance evolution expected in OCDM models, can therefore 
not be ruled out. A similar result is obtained by Barriga et 
al. (2001) who consider primordial step-like features in the 
context of phase transitions during inflation, however their 
features are on larger scales in the linear regime and thus 
would also leave traces in current epoch measurements of 
the CMB and large scale structure. 



6 CONCLUSIONS 

In this paper we presented a series of simulations all based 
on the COBE normalised SCDM model but with Gaussian 
features added to an otherwise unperturbed power spectrum 
(called 'feature' models). Such bumps (or dips) might natu- 
rally arise from non-standard inflationary theories and our 
main purpose was to investigate their influence on the large- 
scale structure of the Universe as measured via galaxy clus- 
ters. We analysed the evolution of the dark matter power 
spectrum, the large-scale velocity field represented by the 
velocity distribution of galaxy clusters, the evolution of the 
cluster abundance, the halo-halo correlation function, and 
density profiles of clusters. We furthermore compared our 
numerical results to analytical predictions based on the 
Press-Schechter formalism (Press & Schechter 1974). 

When comparing the modified SCDM models to a fidu- 
cial OCDM model, we can see that when choosing an ap- 
propriate scale for the added feature the histories of clus- 
ter abundance evolution might be indistinguishable. How- 
ever, there is still a discrepancy in the overall normalisation 
left at the present epoch, with only the slopes of the clus- 
ter abundance evolution coinciding. But as observations of 
the most distant massive clusters only constrain the cur- 
rent abundance to within one or two orders of magnitude 
(cf. Bahcall & Fan 1998), an SCDM model including such a 
broad feature cannot be ruled out. 

Finally we remark that whatever the origin or nature 
is of these bumps and dips, their effect on the evolution- 
ary history of the SCDM model is much more moderate 
than changes in the cosmological parameters of the model 
itself, e.g. lowering the normalisation from as = 1.18 (COBE 
normalisation) to as — 0.52 (cluster normalisation). It is 
therefore an observational challenge to find traces of such 
features, and the best place to search for them might be 
the cluster-cluster correlation function which appears to be 
highly biased with respect to the 'normal' SCDM model. 
This is due to an unusual ratio of power for high- and low-fe 
waves, which will be reflected in the aforementioned 'more 
than biased' cluster formation scenarios. 
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